Multi-particle excitations and spectral densities in quantum spin-systems 

Christian Knetter, Kai P. Schmidt and Gotz S. Uhrig 
Institut fiir Theoretische Physik, Universitdt zu Koln, Ziilpicher Str. 77, D-50937 Koln, Germany 

(February 1, 2008) 



O 

o 



in 



CD 

i 



I 

o 



> 

in 

o 



cd 



o 
o 



X 



The excitation spectrum of the 2-leg S=l/2 Heisenberg lad- 
der is examined perturbatively. Using an optimally chosen 
continuous unitary transformation we expand the Hamilto- 
nian and the Raman operator about the limit of isolated rungs 
leading to high order series expansions allowing to calculate 
spectral densities quantitatively. The 2-particle sector is ex- 
amined for total momentum fc = 0. We show that triplet- 
triplet interaction gives rise to a band splitting. 

PACS numbers: 75.40. Gb, 75.50.Ee, 75.10.Jm 



Antiferromagnetic quantum spin systems are investi- 
gated intensively both in theory and experiment. At 
T = 0, the so called spin liquids exhibit a non magnetic 
short ranged ground state and gapped excitations with 
unusual dynamics. The spin liquid phase is favoured by 
frustration, low dimensionality and low spin. 

Perturbation theory directly allows to obtain informa- 
tion in the thermodynamic limit and is not afflicted by 
finite size effects. The perturbative Continuous Unitary 
Transfomation Method (CUT) introduced previously |^] 
proved to be a versatile tool to calculate high order se- 
ries expansions of low lying excitation energies in spin 
systems (see e.g. [Ij-Q). In this article, we present quan- 
titative results for the spectral densities of multi-particle 
excitations. 

The system considered is the 2-leg S=l/2 Heisenberg 
ladder of which the Hamiltonian reads 

H{Jl, J2)=y ^ [Jl^l,i^2,i +J2 (Sl,iSi_i+i + S2,iS2,i+l)] , 
i 

(1) 

where i denotes the rung (coupling Ji) and 1,2 the leg 
(coupling J2). This system is realized in a number of 
substances M one of them being the superconducting 
compound Sro.4Ca13.6Cu24O41.84 §]■ 

For a given observable O the T = spectral density is 
given by 

(2) 

■00 is the system's ground state. Using the limit of iso- 
lated rungs X = J2I J\ = 0, we map the full Hamilto- 
nian H{Ji, J2) continuously and unitarily onto an effec- 
tive Hamiltonian Hcs{x), x > 0, which conserves the 
number of triplets on the rungs [iJcff(a;), i/(Ji, 0)] = 0. 



The effective Hamiltonian Hf,g{x) is given as an exact 
operator series in x truncated after order 14 in the 0- 
and 1-triplet and after order 13 in the 2-triplet sector. 
For more details see Ref. |l|. The same transformation 
used to derive HcS is applied to observables resulting in 
effective observables Ocs, also given as an operator series 
truncated after 10*'^ order acting on |0), the ground state 
of the effective model (details in Ref. |0]). 

With the triplet vacuum |0) the ground state energy 
Eq is a series in x, Eo{x) = (0|ifcff(a;)|0). Similarly 
the l-(2-)triplet energies are obtained from the states 
\i) {\i,j)) denoting the eigenstates of H{Ji,0) with one 
(two) triplet(s) on rung i (and j) [Q. Result for the ener- 
gies have also been obtained by cluster techniques [P,pl. 
To calculate the Greens function in Eq. (||) for the ef- 
fective operators we utilize the continued fraction tech- 
nique |10| with proper termination 1 1 1 1 . In combination 



with standard series approximations our method allows 
an accuracy within 5% for a; w 1 
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FIG. 1. The 2-triplet spectral densities at total momen- 
tum fc = and X = 1 computed from Eq. (2) with O — Or 
(black) and a test operator (gray) with triplet-triplet inter- 
action switched on (solid) and off (dashed). All curves are 
scaled to give total weight 1. (Further details see text.) In- 
set: 1-triplet dispersion. 

In Ref. iQ we show the fc-resolved spectra for various 
observables at x = 1 in great detail. Relevant for Raman 
scattering [^ 14| is the S = operator Or = Si^iS2,i, 
locally exciting two triplets from the ground state |0) on 
rung z, at total momentum fc = depicted as solid black 



line in Fig. 1 for x — 1. Particularly surprising is the ap- 
pearance of a mid band singularity (MBS) at an energy 
of about 3.4 Ji splitting the 2-triplet continuum into two 
bands. The CUT method allows to identify the parts of 
the 2-triplet energies that solely arise from 1-triplct ki- 
netics 0]. The remaining part is the pure triplet-triplet 
interaction energy. The dashed black and gray lines in 
Fig. 1 (denoted by "no I") depict the spectral density 
of Or at fc = with the interaction part set to zero. 
In both cases the MBS disappears, whereas it appears 
when the interaction is switched on (solid lines "with 
I"). Hence the band splitting is due to the 2-triplet in- 
teraction. To rule out a possible matrix element effect 
we calculate Eq. (||) for an artificial S — Q test operator 
at total momentum fc = which creates the two triplets 
on adjacent sites. The gray ("naive") lines in Fig. 1 show 
the results with / without interaction (solid / dashed). 
Comparing with the corresponding black lines (results 
for the "true" operator Or) one finds that Or merely 
shifts weight towards the lower band edge, being clearly 
not responsible for the MBS. The inset of Fig. 1 shows 
the 1-triplet dispersion Loik) |12] which has a shallow dip 
at small fc-values. This leads to a high density of states 
in the 2-triplet sector at an energy of 2 x w(0) = 3.6Ji. 
Here the dashed gray line displays a van-Hove singular- 
ity reduced to a finite value due to the remaining intrin- 
sic hard-core interaction: two triplets cannot occupy the 
same rung. Due to this large DOS the system becomes 
extremely sensitive to an interaction induced orthogo- 
nality between O|0) and eigenstates \uj) with energy k, 
3.6 Ji. Indeed, the MBS remains detectable even if the 
interaction is lowered to 1% of its full strength. In sum- 
mary we note that the CUT method allows to compute 
spectral densities of quantum multi-particle excitations 
in great detail. The nature of the structures revealed can 
be explained giving insight in the complex dynamics of 
the spin ladder system. 
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